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Abstract 

Groundwater salinity level is one of the most important considerations for 
monitoring soil degradation that threat some regions in the southwest of Basra 
Province, Iraq. The aim of this research is to assign the appropriate and 
effective image processing techniques to be implemented for monitoring, and 
then to evaluate groundwater salinity level map. Landsat TM 2000 and ETM 
2015 images respectively have been selected, as well as ancillary data of the 
available salinity field measurements have been used. Spatial overlay analysis 
between salt affected areas and water table were made to assess spatial 
distribution as well as relationships with these features. The result shows that 
about 47.8% of the areas were low-saline in 2000. This gradually decreased to 
35.9% in 2015. Large area change was observed in the slightly and moderately 
saline soil categories. High sensitive areas of saline soil were about 13.3% in 
2000 and this increased to 19.4% in 2015. Areas highly vulnerable to 
salinization were related to the groundwater salinity level that normally 
occurred on the soil sediment in this location of study. 
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1. Introduction 

Groundwater salinity level is one of the main problems of soil degradation; it is an 
environmental hazard that causes losing the agricultural productivity. Surface 
salinity processes are highly dynamic, so using multi-date images are a suitable 
way to detect the changing state of soil, as well as the technology development 
and extent of environmental change over the last 20th century has given a new 
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Nomenclatures 

E Longitude, degree 

N Latitude, degree 

pH Potential of hydrogen 

TH Total hardness, mg/L 

TDS Total dissolved solids, mg/L 

SO 4 Sulfate or sulphate, mg/L 

Na Sodium, mg/L 

Ca Calcium, mg/L 

NO 3 Nitrate, mg/L 

Abbreviations 

ETM Enhanced Thematic Mapper 

LAO Lood and Agriculture Organization 

GIS Geographic Information System 

MLC Maximum Likelihood Classifier 

RS Remote Sensing 

SI Salinity Index 

TM Thematic Mapper 

WRS Worldwide Reference System 

urgency for monitoring this change. Wide ranges of processing techniques are 
available to discriminate the specttal response in regard to the different soil 
quality [1-3]. 

There are numerous studies on groundwater-associated salinity, but more 
information is required on the effects of groundwater dynamics on soil 
salinization. Lrequent irrigation with small quantities of water is effective to 
reduce soil surface salt accumulation induced by saline shallow groundwater. Salt 
excess in soils has detrimental effect on crop yields and agricultural production 
due to poor land and water management, and results in substantial losses of arable 
soils, especially in the arid and semi-arid areas [4], Lurthermore, salinity also 
affects other major soil degradation phenomena such as soil dispersion, increased 
soil erosion, and engineering problems [5]. 

However, soil salinity is quite time and space dynamic as salinization is the 
consequence of different complex processes of salt redistribution that depends on 
natural conditions, system features, agricultural practices and drainage 
management. In addition, observing the returns and benefits of drainage takes a 
long time (often more than 25 years) that instantaneous measurements of salinity 
do not reflect current conditions [ 6 ]. Observations of many irrigated areas in the 
world also shows that water logging and salinization typically appear only 10-50 
years after the beginning of the project, depending on the initial depth and 
recharge rate of the water table and on drainage conditions. 

Large scale and multi temporal studies of salinity, especially long term 
changes in salinity help to understand the nature of salinization and to evaluate 
the effectiveness of salinity conUol practices. Remote sensing and GIS techniques 
have become tools for the purpose of identifying and classifying saline soils [3, 
7]. This technique is efficient, cost effective, fast, labor saving and accurate for 
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delineating salt affected soils. The integration of remote sensing data, in the form 
of satellite imagery, with GIS has boosted up the ability of delineating and 
mapping soil salinity. Several studies have proved that remote sensing is a 
promising method to identify salt affected soils, especially those with moderate 
and high salinity levels [8, 9]. 


845 Materials and Methods 
847 Study area 

The study area located in the southern parts of Iraq, lies within longitude 
47° 43' to 47° 45' E and from latitude 30° 18' to 30° 21' N with the total area: 
17.792 sq.km (Fig. 1). The soil of Iraq is considered as sedimentary soil, 
especially in the central and southern parts. The annual humidity is less than 50% 
and remains less than 30% during the daytime. The average evaporation exceeds 
2450 mm/year with average annual rainfall less than 100 mm. 


46°30 , 0"E 47°0O’E 47°30’0*E 48°0 , 0"E 48'30*0“ E 



Fig. 1. Location of study area in the southern part of Iraq. 

2.2. Methods 

This research used different analyses, which involves remote sensing for salinity 
change detection, GIS assisted spatial modeling, regression analysis and finally 
validation and comparison of the methods. It involves the integration of thematic 
layers such as geological information, elevation, groundwater salinity level (Table 
1), soil texture and vegetation density in mapping soil salinity. 
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Table 1. Some physical and chemical properties of groundwater. 


Well 

No. 

E 

N 

EC 

TDS 

Cl 

SO 4 

Na 

Mg 

Ca 

NO 3 

1 

47.74401 

30.31809 

11300 

8528 

2850 

2000 

2510 

230 

592 

51.19 

2 

47.74331 

30.32003 

7900 

6000 

1805 

1500 

1550 

163 

544 

41.29 

3 

47.73203 

30.33773 

7580 

5926 

1710 

1100 

1500 

528 

800 

22.2 

4 

47.73604 

30.33841 

15750 

13000 

4750 

1500 

4500 

579 

1280 

51.49 

5 

47.74075 

30.32178 

8510 

6688 

1353 

1400 

1700 

235 

562 

53.79 

6 

47.74677 

30.32113 

9210 

7226 

2233 

1500 

2130 

221 

608 

29.3 

7 

47.74558 

30.31972 

8770 

5948 

1552 

1800 

1400 

197 

530 

16.63 

8 

47.74912 

30.30631 

7190 

5902 

1455 

1600 

1200 

207 

562 

13.48 

9 

47.72943 

30.33897 

7940 

5610 

1710 

1150 

1500 

207 

546 

51.08 

10 

47.74113 

30.32837 

9420 

6628 

2090 

1450 

1800 

160 

499 

29.95 


Landsat images of TM 2000 and ETM 2015 were acquired for the study area 
(Table 2). Layer staking was made for all six bands excluding the thermal band. 
Geometric correction for the three Landsat images was accomplished using a 
topographic map of the study area [9]. Landsat images were registered to the 
topographic map using control points, which were easily recognizable on the 
satellite image. Image enhancement was made to improve the interpretability of 
the images. 


Table 2. Remote sensing data that used in the study. 


Location 

Satellite Sensor 

*WRS 

Path/Row 

Date 

Al Basra- Al- 

Landsat-5 TM 

166/40 

July 10,2000 

Zubair 

Landsat-7 ETM 

166/40 

July 18,2015 


848 WRS is the Worldwide Reference System (WRS), which is a global notation system 
for Landsat data. 


In this research, a maximum likelihood classifier (MLC) was used to retrieve 
urban cover areas water bodies, bare land, vegetation density, and sand land area. 
Because the assemblages are mixed in the class to be mapped, they will generally 
be referred to as land use/cover class. Supervised methods proved superior and 
are covered below. Density maps were created by correlating several test land 
use/cover detection to ground data from published classification map about the 
location area. Urbanization properties, water bodies, bare land, vegetation density, 
and sand land area can be delineated on an image to gather the spectral response 
for similar areas in the rest of the image. By gaining a prior knowledge of an area 
on the image to be classified, response across all bands was matched using 
supervised methods to produce desirable output classes. Salt affected areas were 
clearly identified from other features by higher reflectance in many bands. 
However, it was difficult to differentiate salt affected area from sand soils. 
Salinity Index (S.I.) proposed by Tripathi et al. [10] was applied, which gives 
relatively good results in the re-classification of salt-affected soils. The S.I is 
calculated as 

S.I = (Band] x Band 3) 72 (1) 

where, bandl and band3 represent the spectral bands of the Landsat images. 

Therefore, sample sites identified as salt affected were used for image 
analysis, followed by site verification. Topographic maps of 1:50,000, were used 
in digitizing thematic layers for the overlay analysis. A soil map of 1:10,000 and 
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geology map of 1:250,000 of the area were obtained from FAO soil classification 
and Geological Survey of Basra Province, respectively. To assess the spatial 
distribution of salt affected area with respect to groundwater table. For the 
purpose of interpolation, 10 wells and soil sample data were selected and 
interpolated to generate a continuous surface distributed throughout the study area 
were used. A total of 10-point data with coordinate references were acquired and 
converted in to point shape file in GIS environment. For the spatial overlay 
analysis, the groundwater table point data were interpolated to generate 
continuous surface water table using spatial analysis tools in GIS environment. 
During interpolation, the raster data were sampled to 28.5 sq.m cell size to be 
made compatible with other layers in the analysis [11, 12]. After interpolation, the 
continuous surface water table data were reclassified into three classes based on 
its contribution to the salinization process [1, 4, 13]. According to the soil salinity 
data, extents of soil salinization were graded as series of > 5% (I), 2%-5% (II), 
0.5-2% (III), < 0.5% (IV), which denote high salinization, medium salinization, 
low salinization, unaffected, respectively (Table 3) [11-13]. All the thematic 
layers were generated in GIS environment at a scale of 1:250,000. The software’s 
packages used for this study were (ERDAS ver. 9.1) and GIS (ArcGIS ver. 9.2). 


Table 3. Land surface features of pure pixels. 


Surface 

features 

Biomass 
(kg/30 m z 

Soil moistur 

(%) 

Soil salinit 

(%) 

Salinity 

grades 

Sand area 

0.09 

0.1 

0.05 

IV 

Sparse vegetation area 

14.5 

8.9 

1.96 

III 

Luxuriant vegetation area 

26.6 

18.5 

0.42 

IV 

Wet salty crust 

1.3 

21.9 

9.18 

I 

Dry puffy salty crust 

2.2 

4.6 

7.94 

I 

Salty meadow 

4.9 

5.8 

3.67 

II 


Many models and indicators exist that are used to analyse the magnitude, rate 
and trend of land degradation risk [14-16]. In order to quantify land use/cover 
changes and their regional disparity. Land Degradation Risk (LDR) was 
calculated using the following equation [17]: 

n 

LDR = 2X 

y 

where LDR|. stands for annual detection of land degradation for study area from 
time tj to t 2 , S; is the area of land use/cover type at time tj, Sj.j is the change in 
area for a land parcel whose cover type changed from i to j between 11 and t 2 (I,j= 
1,2,3,4 i4j); n represents the total number of land use/cover types which is four. 
Able to show dynamics of land use/cover changes and their variation across 
different regions [18], 

The land use/cover change map derived from the satellite image was digitized 
and edited in ArcGIS. Other derived attribute data, viz., and population pressure, 
were incorporated into the GIS database using FoxPro . Final tally of land 
use/cover change for study area was determined by merging all ranks of change 
associated with the land use/cover indicators. Prior to the merging, each rank level 
was converted into a numerical value according to an established linear and 
continuous mathematical equation. The area of land use/cover detection at each 


Zs; X (t 2 — 6 ) 1 xl00% 


( 2 ) 
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overall level of change within study area was ascertained by overlaying the final 
land use/cover change map with the administrative boundary map using ArcGIS. 


3. Results and Discussion 

Table 4 shows significant changes in the LULC can be recognized and type of 
LULC conversion taking place can be identified. From 2000 to 2015 sand land, 
urban areas, and bare lands saw a relatively dramatic increase. The areas that 
contributed the most to this change were vegetation land. This may suggest 
logging and development. During the same time span, urban areas also increased 
in size. The majority of this change came from the development of vegetation 
land into an urban class. Water bodies saw a decrease 10.6% in size. Some of the 
water bodies converted to unused land, while some of it was converted to urban. 
Vegetation areas also saw a decrease in size. Most of it was converted to sand 
cover area, while a smaller but significant portion was developed into urban areas. 

Table 4. (LULC) classes monitored from satellite image for the study area. 


(LULC) classes 

Area 2000 

Area 2015 

Amount Change 

Vegetation land 

km 2 

% 

km 2 

% 

km 2 

% 

Sand land 

5.32 

29.9 

4.68 

26.3 

-0.64 

- 12.0 

Urban area 

5.64 

31.7 

6.03 

33.9 

0.39 

6.9 

Unused land 

1.83 

10.3 

2.12 

11.9 

0.29 

15.8 

Water bodies 

2.75 

15.5 

2.95 

16.6 

0.20 

7.2 

Vegetation land 

2.25 

12.6 

2.01 

11.3 

-0.24 

- 10.6 


Two scenes for each of Landsat TM and ETM+ for the period 2000 and 2015 
respectively have been used to create an image index of Salinity Index (SI) 
depending on the spectral response of the soil. This index reveals an increase in 
salt-affected area during the mentioned period. Ancillary data is accurate strategy 
that indicators different salinity level in the term of spectral response. The results 
show the incorporation of the EC relationship with Na + , Ca ++ , and Mg + content 
which provides a better and more objective indicator of soil salinity and indicates 
that there is greater soil salinity with increasing Cl, SO 4 , and NO 3 content. 
Accordingly, many information layers of the salinity elements have been created 
as shown in Figs. 2 and 3. As well as there is good contribution between EC and 
spectral response in Landsat images specifically in Band 1 and Band 3 which is 
considered the most accurate band to detect and evaluate soil salinity. 

The salinity index (SI) was used to enhance the saline zones and suppressing the 
vegetation, which was calculated as the ratio of the difference of red band to NIR 
band. The values of the salinity index (SI) ranged between 1 and -1, from these areas 
with value 0 and less than 0 were classified as none saline. About 47.8% of the study 

area was low sensitive areas in 2000, which gradually decreased to 35.9% in 2015 

2 

(Table 5). Out of this, 8.50 km - under low sensitive areas in 2000, whereas the rest 
area was changed in to slightly saline. Moderately saline soil was about 38.9% in 
2000, which was increased to 44.7% in 2015, respectively. 
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Fig. 2. Information layers of Na, Ca, and Mg as contours line. 


Table 5. Salt affected area class derived 
from SI and change rate for 2000 and 2015. 


Salt affected area class 

Area 2000 

kin % 

Area 2015 
km z % 

Rate of Change 
km 2 % 

Low sensitive areas 

8.50 

47.8 

6.39 

35.9 

-2.11 

-11.9 

Moderately sensitive areas 

6.92 

38.9 

7.95 

44.7 

5.80 

32.6 

High sensitive areas 

2.37 

13.3 

3.45 

19.4 

6.10 

34.3 

Total 

17.792 

100 

17.792 

100 

- 

- 


The results indicate that salt affected areas were distributed throughout the study 
area. In the central and in many of the northern parts of the farm was moderately 
saline. As illustrated in Fig. 4, the north-western parts of the location were of high 
salinity. A map of salt affected soil for the two different study years revealed that 
moderately and slightly saline soils were concentrated in the central area of the study 
location. Figure 4 shows the area ratios and statistics characteristics of soil salinity of 
different groundwater depth types in the study area. There are clear 
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differences of soil salinity among them. As soil salinity problem was closely 
related with salt water concentration, the salt affected areas mapped were overlaid 
with the groundwater table depth. The Groundwater table depth of the area was 
within 10-30 m below the ground surface. About 54% of the salt affected area of 
groundwater salinity level. The classification of groundwater table depth with 
different level class and risk salt water concentration indicated that groundwater 
salt content is the major problem of the occurrences of soil salinity especially in 
critical areas. From the prediction model, the total area identified as salt affected 
soil was 25%. Out of this, 11% is lying within critical areas and 64% was lying in 
potentially area classes of the total salt affected soil exists (Fig . 4). 



Fig. 3. Information layers of Cl, SO 4 , and NO 3 as contours line. 


Since salinity is a dynamic process it is important to monitor salinity process 
and map its spatial distribution regularly. Although geostatistical technique is 
available to see the spatial structure it takes lots of effort in collecting sufficient 
soil samples and their laboratory analysis (Fig. 5). In this respect combination of 
geopedologic interpretation of the area and incorporation of band rotation seems 
to be useful in accessing salinity problem quickly. Results, using spatial analysis 
methods, showed that 11.2% of land had no risk of land degradation by soil 
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salinity, 26.1% had moderate, and 62.7% of the total land area was at a high risk 
of land degradation by soil salinity. In conclusion, the study area is exposed to a 
high risk of soil salinity. The results of the statistical analysis showed that saline 
area has a significant correlation with vegetation cover negative change (0.89). 



Legend 

Salt affacted area 
Critically area 
Potentially area 


□ 

n 
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N 

A 
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012 4 KM 


Fig. 4. Created salinity index (SI) for the periods 2000 and 2015. 


The area of salt affected soils that resides within 100 m and 300 m distance of 
wells are found to be only 40.8% and 29.3%, respectively (Fig. 5). The result 
indicated that poor irrigation and water management systems and small canal 
structures lead to secondary salinization. Secondary salinity resulting from 
modern irrigation occurs due to accelerated redistribution of salt in the profile by 
high water table. Irrigation was the main cause of rise in shallow groundwater 
table under intense evapo-transpiration conditions and this led indirectly to soil 
salinization. In general, stretching from the wells buffer, soil salinity shows 
decreasing tendency with the decrease in distance from the wells buffer. It reflects 
the significant influence on soil salinity by the salt water around. . The statistical 
analysis showed this index (SI) has a significant correlation with water bodies’ 
positive change (0.92). 
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Fig. 5. Geostatistical distribution of salt affected area in the study area. 

Factors causing soil salinity include inappropriate and excessive irrigation 
without adequate drainage system, irrigation water quality, rising saline water 
table, climate, rainfall history, local topography and soil composition and farming 
practices. Therefore, increasing soil salinity at the surface is most likely to vary 
according to the distribution of these different factors across the landscape. It was 
revealed that the spatial distribution was not highly influenced by the features 
considered except groundwater table. The results demonstrate that modeling and 
mapping spatial variation of soil salinity based on remote sensing data is a 
promising approach. Based on this view, to assess the distribution of saline soils, 
the water table map was overlaid with a map of salt affected soils. The rising 
water table brings salt from deep in the soil up to the surface, causing salt 
accumulation. The present results also suggested that a rising water table and salt 
accumulation at the surface combined with a high evaporation rate are likely 
factors that have resulted in the spatial variation in soil salinity of the area. This 
study shows how regression analysis, coupled with remote sensing images, could 
successfully predict and map spatial variation in soil salinity over an area 
vegetated. Thus, the information presented can help farmers, scientists and 
engineers to manage soil salinity problems affecting the ecosystem. Additionally, 
the simplicity of this approach, with its satisfactory accuracy, can contribute to 
soil salinity prediction and mapping at lower costs than conventional approaches. 
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Figure 6 and Table 6 shows the distribution of land degradation risk areas in 
study location, it is clear that the sensitive areas to land degradation in study area, 
where the soil quality, climatic quality and management quality are low; these 

2 

areas represent 12.8 % ol the study location area (i.e. 2.28 km ). The areas of high 

to moderate sensitive for land degradation exhibit the rest of the parts of the study 

2 

location as it represent 87.2 % of the total area (i.e. 15.51 km ). The northern 
parts of the study location are characterized by a high risk for land degradation as 
they represent 61.9 % of the total area (i.e.11.01). The low risk for land 
degradation is due to the good vegetation cover and soil quality. The results of 
this study indicated that land degradation results from natural and anthropogenic 
factors. Overlay of land degradation processes layers interpreted from multi¬ 
temporal remotely sensed materials in a GIS, in conjunction with field 
investigation, revealed that the spatial extent of sandy desertified land in the area 
has drastically expanded during the fifteen-year study period (2000-2015). 



Fig. 6. Land degradation risk assessment in study area. 


Land degradation processes in the study area was assessed through 
consideration of both natural (vegetative index, soil index, climatic index, drifting 
sand) and anthropogenic (Land use change) factors in the study. It was found that 
most of the study locations were highly land degradation. The overall sensitivity 
of land degradation change has worsened during the study period with degraded 
areas accounting for 61.9% of the total area in 2015. There is a clear trend in the 
spatial distribution of the land degradation direction within the study area, which 
goes from the Northwest toward the Southeast. The risk has risen considerably, on 
an average, by 40% for all western parts of study location between 2000 and 
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2015. In particular, the risk has increased considerably for those areas not 
previously considered highly vulnerable to degradation. Consequently, the 
disparity of land degradation hazard among the study locations has shrunk as all 
of them are at a higher risk in 2015 than ever before. The accentuation of land 
degradation is attributed to conflicts among human interest, limited land resource, 
and fragile ecosystems. Inappropriate human activities such as excessive 
exploitation of natural resource and mismanagement of land, to a certain extent, 
have contributed to the land destruction. 


Table 6. The categories of land degradation 
and the proportion of each category 


Class 

Area (km' 1 ) 

% 

Low land degradation areas 

2.28 

12.8 

Moderately land degradation areas 

4.50 

25.3 

High land degradation areas 

11.01 

61.9 

Total 

17.792 

100 


4. Conclusions 

This study proposed a system depends on satellite images as the most important input 
data to the GIS, and it is ready to be imposed by other promoted required data for 
updating or any other environmental study in order to reduce the expense that is 
needed to construct other related study. It provides a tool more easily to manage the 
immediate attention with greater monitoring. Remote sensing and GIS based 
integration and spatial overlay soil salinity model using multivariate analysis results 
showed that most of the moderately saline soils are on saline water table depth. 
Intensive irrigation practices and excessive use of water and poor irrigation 
management can cause soil salinization. Areas where groundwater table lay above 
critical depth increased as a result of higher evaporation and this has contributed 
directly to the occurrence of moderate and slight salt accumulation in the area. The 
model has good predictability as it considers most of the factors, which contribute to 
soil salinization. Modeling and mapping spatial variation in soil salinity based on GIS 
analysis and remote sensing data are promising approaches, as it facilitates timely 
detection with a low-cost procedure and allow decision makers to decide on necessary 
actions, to be taken in the early stages to control soil salinity and to conserve 
agricultural lands and natural ecosystems. 
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Abstract 

The insertion of the pre-reinforced element allows to overcome the 
initial stresses in the reinforced structure, reinforced under the operational load. 
In our time, a rich practice was gathered in the area of preliminary 
intensification of bending and tightening of elements. Along with this, there are 
only a few studies on the preliminary intensification of such tightened elements. 
Preliminary intensification of molded building elements can be carried out in 
various ways (with the help of pre-strengthened tugs, telescopic pipes, 
intelligent communications, etc.). The reinforcement of the rods pressed in 
metal structures with the help of pre-strengthened tugboats allows to reduce the 
volume of construction and installation works to a minimum. 

Key words: initial tensions, reinforcing, reinforced, construction, 
stability, pre-intensified, diaphragm, moment, 
resistance, elastic supports, collar, matrix, power 
crisis, protection coefficient. 


Introduction 

The investigations showed that the area of the cross-section calculated taking into 
consideration combined work of the reinforcing and reinforced elements is 
significantly less than the one calculated taking into account separate works of the 
above mentioned elements. 

In the reinforced elements, initial tensions are to be taken into the 
consideration in calculations. Reinforcement process in the constructions under 
the loads should comply with Construction Norms and Rules of on the basis of 


N r 0,6N 0 


condition. Here. 


moment of reinforcement; 


N, 


is a force in the previous element in the 
0 - is an account force in the pivot before 


reinforcement. It is due to the fact that while 1 > ’ 0 heating of the pivot (at 

the result of welding) leads to increase of bending and loss of stability in the 
pivot. In order to increase efficiency of the reinforcement under the load, pre¬ 
intensification of the reinforcing elements are realized. Pre-intensification of the 
reinforcing elements allows to increase elastic work limit of it. [7,8]. 

Reinforcement of load-bearing constructions and the ways of regulation of 
tension in them. Pre-intensification of the reinforcing elements can be realized by 
means of electrothermal and electromechanical methods. Besides it, hydraulic 
jacks and other force using methods can be applied for pre-intensification. The 
reinforcing element is pressed (intensified) in advance with the help of any 
method, and then it is connected with the reinforced construction reliably in the 
ends. The reinforcing element released from the tows will begin to be tightened 
and so, will reduce the tension of the reinforced construction (tightened pivot, 
pillar) to some extent. The reinforced pivot is transversely connected with the tow 
by means of diaphragms and it provides their joint effective work. 
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The gap between the diaphragm and the tow (very often met in practice) 
affects the stability of the pivot. Account scheme of such a tow can be considered 
the pivot with elastic intermediate support. In the middle pass of a one-diaphragm 
element, in intensification process a curve emerging at the moment of lock of the 
gap between the diaphragm and the pivot can be determined in the following way 
[ 2 ]. 


/i=/o+^=— ( 


N, 


750 N-N, 


■) 


Here, 


/o ; 


is an initial curve; 


Sy. 


( 1 ) 

is a gap between the tow and the diaphragm. 

We 


fl = 


On another hand, if we suppose 
determine that: 


(N e -N?) 


from [2], we can 


/t - /o + - 


f 0 N e 


(N e -N?) 


( 2 ) 


N° = 


M 

fo +&1 


Here, 

(3) 

N° 

1 - is a longitudinal force emerging at the moment of the lock of the gap 
between the diaphragm and the tow. 

If to approximate the tilted axes of the pivot with sinusoid 
7ZX 

y = fsm (—) 

, we can write the following for the locking moment of the gap 
(between the diaphragm and the tow): 


7tx 


Ely" = —N j/j sin — 


of 


Here, we can write the following, if we take into account the values 
from (3): 

Ely" = -NJ sin — = -N e ^ f x sin — = -NJ l sin — 

* /i * * 

(4) 

Increase of the intensified pivot will stop after its distortion; the reason for it 
is an elastic reaction of the tow. 

Elastic reaction of the tow can be determined from the equilibrium condition 

both in elastic and plastic stages. Under the impact of the d force, the tow will 
distort from its rectilinear position and a curve will emerge in the middle pass ( 


^ 1 ). Such a curve emerges at the result of increase of ^ force up to ^ 


force. 


According to the equality condition at the place of dislocation of the middle cuts 

f 

of the tow and the pivot, the curve of the tow will be equal to J 2 (Figure 1). 
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Fig. 1. Account scheme of the pivot pre-intensified by means of the longitudinal 
tow 


f = I AO 2 — AO 2 

From AOOl triangle (Figure 2) we can write J 2 » 1 

AO = 0,51 d + M d . AO, =0,5£ d 

So, we get: 


f 2 =0,5l(£ d + 


(5) 


= o 5f 


N 


( E d A.) 


„ N 2 

+ 2 - 2 


E d A d 


If we take into consideration approximation (at the result of shortening and 
bending of the pivot) of the ends of the intensified elements [3,4]. 


f 2 =0,5(1 d -S m -A m \ 


N 2 


(E 2 A 2 ) E d A, 


- + 2- N2 


d ls -d 


( 6 ) 


Here, m - is the value identifying the shortening of the pivot. 

N 

w force will increase as much as 2 
reaction emerging in the tow in the middle pass. 

'e are uni 

f 2 (N 2 ) 


N R 

The tow force will increase as much as 2 at the result of the impact of d 


N f 

As 2 force and J 2 curve are unknown, equation (6) should be solved 


together with another dependence 

R d =Qd 

N „ 


Sy 


. Also, elastic resistance of the tow is 

r =n 

unknown 

Dependence for the added '' 2 force after gap is closed can be written 

from equilibrium condition of internal and external forces (Figure 1). As the @ 

angle is small ( COS @ ^), this equilibrium condition can be written as follows 

after some transformations: 

N -N 

N,f,cosa - RW ■+—-W = - -W - N 2 f,cosa 

A A 


N 2 =- 


Here we get: 


N, (/[ cos a + p)~ RW 
f, cos a + p 


(7) 


R W 

Here, ’ - are account resistance of the metal and resistance moment ol the 

cross-section cut of the intensified element. 
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cos a = 


E d A d 


N 2 + A d E d 


( 8 ) 


If we write this formula in the formula (7), we will get: 

Aye/, + p)-RW _ Kl/j +p)~ rw\n 2 + E d A d ) 


N 2 < 


(9) 


f x E d A d /(N 2 + E d A d ) + p f\E d A d + p(N 2 + E d A d ) 


When the sum of ^ d , and ^ is known, it is possible to identify the value 
ike 

fi 


N N 

of 2 . If we take into consideration the value of 2 in (5) and (6) formulas, 


we can calculate 

R — O 

It is possible to determine ^ d force, supposing d i 
the tow as a tightened agile thread from the condition of equality of the moments 
to zero in regard to 0 point (Figure 2). 


Qd fnt-pp cnnnncinrr "d ^ d j n middle paSS Of 


Q d --~N 2 cosa- f 2 =0 


Qd = ^^ 2/2 cos« = 


4 N 2 f 2 


£(l + N 2 /E d A d ) 


Here, 

( 10 ) 

^ — is the distance between the connected ends of the tow. 

It is also possible to determine this value of ^ d from the equilibrium condition 
of the intensified element in elastic stage. 



Fig. 2. Scheme of the reinforcement by means of a tow of a reinforcing element 
having just one diaphragm in the middle pass (account scheme of the tow). 

The moment of the external forces in regard to the point located in the middle 
pass, is regulated with the moment of the internal forces. 


£m 0 = 0,5ft A 


A', /, + AM/, + /,) cos « (A' / A,)]»/ 


From here we get: 

Q d =*[N l f l -N 2 (f l +f 2 )cosa + (R-N/A 0 )W 0 ] = 


4 

7 


^/i-(v 2 (/i + / 2 ) 


E d A d 
+ A d E d 


N 

+ RW -W r 


( 11 ) 

Here, 


n = n x +n 2 . 


R 


- is an account resistance of steel; 
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E M 0 =N 1 f l + N 2 (f x + f 2 ) cos « - ^ = -(.R - )W 0 


From here we can write: 


Q d =- 


N 


NJ , + N, (./, + /, )cos« + (7? - —)W 0 

A) 


A/,+A(/i + / 2 ) 


At A 


N 

-+ (/?-)W 0 

iV, + a, e, a 


under the impact of 1 ? ^2 


'2 1 

Q d 


( 12 ) 


/=/ i +/ 2 


, the general curve J J 1 ' J 2 will be as 
differential equation of the bent axis of the obtained intensifying element: 

1 2 , 


El 


cry 

dx 2 


N 1 f l cosa + N 2 (f l +/ 2 )cos«- 




-cos a 


In this equation, if to write the values of 
will get: 

■ d 2 7TX 

El —y/(sin —) = -(N 1 + N 2 )f l cosa 
dx c 

If to conduct differentiation, we will get: 

Ei f = cos a(N l + N 2 )fi 


y = f sin(^) 


and 


Q d 


(9), we 


i 


(13) 


N„=Ei^ 

e 


If to consider c , we can write: 

N e f = {N X + N 2 )cosa- f { -N^^osa + N 2 f l cosa 

We get the following after some transforms of the equation (13): 


N 2= N e 


(/l+/ 2 -A) 


/t 


A+/ 2 - 


N 

\ — e ' J 2 \ 

E d A d A 


= N e (fS l ) 


E d A d 


f\E d A d —N e f 


(14) 

f N 

If to solve this equation according to J 2 and to write the value of 2 (7) in its 
place and to conduct some transforms, and to suppose the small value of a as 
COSC/ ~ 1 ( we will get the following: 

, [jV,(/, +p)~ W]-/,(l-JV,/E„4,) + ^,/„(/, +p) 

[N,(f 1 +p)-RW}l + N,/E J A J ) 

(15) 

When 2 force is c OS(Z 1, it will be equal to zero and it means that the gap 
between the curved diaphragm determined by the formula (14) and the tow will 
be a curve formed before it closes. 


When A and 
of the formula (10). 


A 


2 are known, it is possible to calculate ^ d with the help 


Q d 
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It is possible to check hardness and stability of the intensified element using 

the values of 1 , 2 , 7 1 , ■' 2 , . According to the account scheme of the 

intensified element (Figure 3), taking into consideration the formula (10), and 

supposing that the force ^ d balances with ^ 1 , we will get (a balance equation 
is written in regard to the middle pass (zero point)): 
cr <R 

y (16) 

cr = (N, +N 1 )cosp\ + M N ^--Q d f --^ ; 


Here, 


A 




4 W 


M n = [NJ\ + N 2 U\ + f 2 )]cos« 


,L 


(-) -/i =(-) z cos^ 

From the formula ^ z - 


we get (Figure 5.3): 


COS /? = ,!- 


4/r 


Taking into consideration the values of & > ^N > COSCC > COS/? an( j ; the 
condition (12) can be written as follows: 

1 


cr = 


N l+ N 2 


i-14i + a/i 


IT U-NJE.A, 


<R 


(17) 



Fig. 3. Account scheme on verification of stability of the intensified element 

if , . , /,// = 1/200 . . cos p«1 

It to take into consideration that even when ^ 1 , it is ' 


w 


N - N { +N 2 


The formula (17) can be written as follows: 

A 


rf 1 + N 2 /E d A d 


(18) 


Or, 

<7 = 1 + ^ 


r'- R ' 


(19) 


r\ l + N 2 /E d A d 
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Here, ^ - is the distance from the center of gravity of the cut until the most 

y 

intensified fiber; 1 - is a radius of an inertia cut. 


As it seems from the formulas (3), (18) and (19), the load capacity of an 
element with one diaphragm in the middle pass intensified by the help of the tow 
will be as much more as n than the similar element without diaphragm (cross 
section) where. 


£ 


n = 


/o + fi 


fo +^i 


N„ 


E d A d 


/,+</„+<5, X1--A-) 

h d A d 


( 20 ) 


Availability of a gap between the diaphragm and the tow leads to reduction of 


n parameter. 

The reinforced element connected with the clamps and other fastening 
means with the reinforced pivot can be considered as a pivot on elastic supports. 

In this case, the reinforcing element undergoes a pressing by a pre-intensifying 

N 

force q , and elastic supports serve as one-sided communication. 

Collapse rate of such a one-sided communication (elastic supports) will 
depend both on bending hardness of the reinforcing (main) element and location 
of these supports alongside the pass. The number of intermediate supports should 
be taken so that besides ensuring stability of the reinforcing element, also the 
necessary level of pre-intensifying force was achieved. 

In other words, the number of the intermediate supports (clamps) should be 
taken so that the bending stiffness of the reinforcing element was minimum [2,3]. 
For the above stated case: 


N„ P 


kr 


( 21 ) 


That is, pre-intensifying force q should not be less than the crisis load kr of 
the intensified element. 

In order to solve this issue, account scheme of the most unsuitable work of 
the reinforced and reinforcing elements should be considered. Fastening clamps 
(elastic supports) should be installed in such a way that they did not limit the 
movement of the reinforcing pivot in longitudinal direction. Exacerbation of the 
reinforcing element in advance can be accepted as a pivot with elastic 
intermediate supports and exposed to longitudinal pressing force (Figure 4). 



Fig. 4. High tide scheme of reinforcing element. 

At the result of the preliminary intensification (longitudinal pressing) loss of 

stability depending on the number of the high tide waves ( ^ ) of the reinforcing 
pivot can take place in the different values of the crisis load. 

When hardness of intermediate supports is small, the reinforcing pivot will 
swell only over one wave. When hardness is increased, high tides of the 
reinforcing pivot will increase and the value of the crisis load will increase. 
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In this case, the number of the high tide waves will be equal to the crisis value and 
increase of hardness of the supports will not lead to increase of their number. It 
can occur in the case when there is enough space between the reinforcing and 
reinforced elements, otherwise, additional high tide waves will emerge in the 
places where clumps have been put. If collapse of the supports is very little, then 

k n 

their protection coefficient “ u ” will be determined as follows: 


■ a- 


7T 2 EI 


( 22 ) 


Here, ^ - the distance between the intermediate supports is given in the work [11] 
of the value of a - coefficient, it depends on the number of the passes and is 
called characteristic numbers of protection coefficients. 

Ensuring stability of protection coefficient (i.e. ensuring stability of the value 

k 

of “ °” coefficient) bears important practical significance as the given value of 
this coefficient being more or less than that one calculated by means of the 
formula (22) indicates availability of zero points on elastic supports during the 
high tide. 

Let’s have a look at an elastic system loaded with external load groups given 
via certain Y parameter. Let’s make small replacements in the system 


8, (i = 1 , 2 


8 


. In this case, besides ° replacements, additional external 
forces (support reactions) will emerge in the considered system: 


R = 


R, 


R„ 


(23) 


Given group of external forces creates moments according to 1 - points: 

^ My — /AS (24) 

Here ^ - is a numerical matrix depending on the type of the device and loading 
regulations. 

On another hand, external reaction is a bending moment formed from R 
forces on 1 - points: 


M r =L m R 


(25) 


Here, m - is a matrix compiled from the ordinates of the penetration lines of the 
moments. 

R - reaction forces depend on 8 . values. 

R-BS (26) 

Here, ^ - is a non-specific matrix. 

So it can be written [9]: 


M,=L m BS 


(27) 


Total bending moment on considered 1 -points: 

M = AS +L.BS 


(28) 
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Let’s replace elastic load ^ M I El e ] ast j c ^ loads ^ 

Then let’s write the equation of curves taking into consideration replacement of 

the bending moments with impact matrix in the girder cuts: 

y = LmW (29) 

Elastic loads for the whole systems are determined with the following formula: 
r QpQm r N..N.. 

GA 


r M M r Q U 

w = + k VM^ds + f- 

" J El J nA J 


EA 


-ds 


(30) 


If to replace the opposite points of the moment diagram between the neighboring 


n — 1 n 


and n cuts with a straight line and to take into the consideration 
only the first integral in formula (30), elastic loads can be calculated with this 
formula: 


W = rDM 

6EI 


(31) 


Here, ^ - is a distance between the cuts; ^ - is a Jacobi matrix; M _j s the 
bending moments in the cuts [10]. 


D = 


a. 


a 


12 


a 


21 


a 


22 


a 


a 


32 


23 


CC 


...a„(n-l) a n 


(32) 


The elements of the matrix are determined in the following way: 

a k (k -1) = p„ . a kk = 2 (p n + p n+l ) . a k (k + 1) = p n+1 

(33) 

Here, ^^ 

(34) 

For the stability issue of the pre-intensified element we have considered: 
Pn ~ Pn+\ ( S n ~ S 0 , I n - 1 q y In that case &k (k ~ 1) - 1 . OC kk - 4 . 

a k (k +1) = 1 

So, we get: 

6 El 


M 


C L 


Here, 


(35) 


• _ ^ ^ 

is the ordinates of the elastic loads on ’ points. 


El 

4 1 

1 4 1 

1 4 1 

.1 4 

Here, -is a modular Jacobi matrix. 


A = 


(36) 
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In this case, we will get the following for determining the replacements: 

y = 


El 


(37) 


If a high-precision is requested in determination of replacements, then in this case, 
it is possible to replace the diagram of the bending moments between two 

neighboring points ( n ~ ^ , n , H +1) with a square parabola passing through 
three points (Figure 5). 


Sn 


Mn+1 


Mn 


1/Sn 


1/Sn 


Mn+1 


1/Sn+1 


1/Sn+1 


Fig. 5. A scheme on calculation of a pressed pivot with intermediate elastic 
supports by the method of elastic loads. 


In this case, more precise value of ' - will be identified in such a way: 

y = 


6 El 


(38) 


( 1 . 0 , 1 ) 


If we will take into consideration that ^ D(A£) 5D (M ) j s j n e q Ua ] 

unit (linear) hardness, then, it can be written [9]. 

6 El (39) 

Here, 

" 1 0,1 


D 


(1, 0,1) 


0,1 1 0,1 
0,1 1 0,1 
.0,1 1 


If to write the value of (28) in its place in the formula (35) : 

™ = WAS + D^BS) 

6EI 


(40) 


(41) 

Using f ro m the impact matrixes of the moments, we can write: 


8 = LW = L„ 


6 El 


(yD t AS + D t L m BS) 
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( 42 ) 

We can write the formula (42) in such a form if to include 

C = LDB 


K = LDA 


, a S/6EI markings: 
8 = ayKS + aCS 


yKS = 


-E-B 

a 


8 


Or , 

(43) 

If to conduct a number of transforms here, we will get: 

(K 0 --E)e = 0 

7 (44) 

Here, K 0 = Kid/a)E - B]-\ 

The analysis of the formula (44) shows that the crisis case of the system will 

emerge when it is & ^ ® . 

In that case. 


K 0 --E 

7 


= 0 


(45) 


This equation becomes possible only if we accept the values of characteristic 

f* 1 / y 

2 >...> n ), for ' sum, the crisis for these 


numbers of \ ^ ^ 


i / n. i / r 2 . - l/ 7„ 


characteristic numbers will correspond to ' 1 ; ' 2 ; . . parameters. 

Finally, solution of the issue leads to determination of the minimum value of the 
crisis load [6]. 

Yi = llR i 


As it seems from the 


y 

equality, the parameter of ' kr corresponds to 


the value of 




. So, for solution of stability the greatest characteristic number 


of 


K n 


matrix should be found. 


We can write the formula (45) in the following form after a number of transforms: 

I K -AE\ = 0 

1 1 (46) 

Here, ^ = 1/a 7, 

From the formula (46) we get: 

6 Ei 

7kr ~ TA 

^^max (47) 


In the formula (46), K matrix is determined as the sum of three matrixes: 

K = L ^ DA (48) 

The reinforcing element works as a pivot resting on an elastic basis and 
exposed to pressing. Here reinforced element serves as the elastic basis. In this 
case, the loads are passed to the reinforced elements by means of clumps and high 
tide semi-waves of the reinforcing element (Figure 6). 


_ 
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Fig. 6. Account scheme of the reinforcing element. 

If to take into consideration that there is a certain gap (at the result of hard 
fastening of clumps) between the reinforcing and reinforced elements and the 
impact of the clumps is like one-sided communication, high tide of the reinforcing 
element will have semi-waves in minimum quantity. We suppose the collapse of 
the intermediate supports located in equal distances in the middle of the 
intensified element equal to each other. And we accept elasticity coefficient 
(hardness degree) of a separate support equal to the reaction force creating a 
single replacement in perpendicular direction of an axis of the pivot of that 
support. General hardness of all intermediate supports is accepted approximately 
as follows: [1,4,5]. 

- K { a { n 

Here, ^ - is the hardness of the most collapsed support; n - is the number of the 

supports; 1 - a coefficient accepted depending on the type of the support 
(jointed or stiff). Impact of elastic supports located in equal distances can be 
replaced by an equivalent impact of elastic environment. 


Hardness coefficient of such an environment can be expressed in the following 
way: 


(49) 


t. 


Here, " - is the distance between the intermediate supports (clumps). 

In this case, the crisis load of the intensified element will be determined by 
means of the following formula [1,5]. 


P kr = 


n 2 Ei 


(Kf + 


K 2 7T 4 Ei 


) 


(50) 


Here, 




is the number of semi-waves formed at the moment of high tide in the 


intensified (intensifying) element; 


El 


and ^ - are consequently hardness and 


length of the element. Supposing hardness of the elastic environment as 


Y 


minimum, we determine such a value of the crisis load that in this value, loss of 
stability of the pivot takes place in two neighboring forms, i.e. loss of stability 

K K — 1 

emerges when the number of semi-waves are equal to 1 and 1 
If to suppose the minimum value of the crisis load as much as the number of the 

A'j-l . 

1 ir 


semi-waves as 


in the formula (50), we will get: 


pK 1 _ El 

^kr ~ n2 




{K.-iy^Ei 


ii 


If loss of the semi-waves is equally possible both in 


then in this case, we get the above shown from the condition 


(51) 

** 1 and - numbers, 

P, K J~ l = P, K 


kr 


kr 


tt 2 EI 


(K l -1) 2 + 


(K 1 - X) n El 


7r 2 EI 


(Kf+- 


K 2 k a EI 


i) 
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Here, after some simplifications we get: 

-4=*, ! (k- if 

n El 

From this formula we identify Y . 
r _Kl(K t - \)7r 4 Ei 


(52) 


1 (53) 

If to write this value in its place in the formula (50), we will get: 


4 = 


Or, 


K 2 El 


( 24 - 24 + 1 ) 


(54) 


(55) 


4 = 424 - 24 + 1 ) 

P 

Here, e - is a crisis force in the girder without intermediate supports. 


If to write the account value of the pre-intensifying force 


4 , P, 


kr 


to the 

left side of the formula (55), we will get the minimal quantity of the semi-waves 
of the intensified element in high tide: 

2 4 

24 2 - 24 +l = ^- 

p 

(56) 

Start of the high tide of the intensified element as a rule, takes place at the result 
of availability of the big spaces between it and the reinforced element (presence 
of an initial bending and not good fastening of clumps). At the result, the number 
of the semi-waves will be minimum and bending points of the intensified element 
will be formed most likely in the cuts where the clumps have been put. 


Q d 


Outcome 

Formulas have been worked out in order to calculate 4 , 4 2 and '^ d 

parameters included to the balance equations and considered as the pivot working 
with pressing and reinforced with pre-intensified element. A rule is given for 

N 

checking hardness and stability ol the intensified element using the values of 1 , 

4 1 fi fi anc j Qd 

It has been noted that the load capacity of the element with one diaphragm located 

fi 

in the middle pass and intensified by means of a tow will be as much as “ ” than 

n = (N, +4)/iV 

the similar element without diaphragm ( 1 z e ). 

A formula has been worked out to determine the replacement at the result of the 
pre-intensification of the tow connected with the main element by means of the 
diaphragms. Besides it, canoenic equation has been compiled in order to ensure 

the stability, and a matrix has been formed for 0 protection coefficient of the 
tow. 

A formula has been worked out for determination of the crisis load of the 
intensified pivot. It was noted that loss of stability takes place over two 
neighboring form and allows to create a formula for calculation of the value of the 


crisis load depending on 


.K , 


’ number of semi-waves. 
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tability of the pivot reinforced by means of the reinforcing elements has been 
considered and it has been noted that the value of the crisis tension of the pivot 
reinforced under the load by means of pre-intensified elements is very close to the 
value of the crisis tension of the unloaded pivot. 
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